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THE KINETICS OF ENZYME INACTIVATION 


ALSTON S. HOUSEHOLDER AND GEORGE GOMORI 
THE UNIVERSITY OF CHICAGO 


The course of hydrolysis of 8-glycerophosphate catalyzed by a group 
of different enzyme extracts, both with and without the addition of Mg, 
with and without preincubation of the enzyme, has been studied and the 
results discussed on the basis of a mathematical analysis. In all the ex- 
tracts, it appears that two distinct and independently acting constitu- 
ent enzymes—or perhaps “principles” of the same enzyme—are present, 
one acting much more rapidly but also more rapidly inactivated than the 
other. Storage in the refrigerator changes markedly the behavior of 
both constituents, though in different ways. There is evidence that in 
some cases an enzyme is limited in its hydrolytic “capacity” in the sense 
that after an enzyme molecule has decomposed a definite number of sub- 
strate molecules, it thereafter becomes entirely passive. Further, there 
is evidence, in the case of one extract, that the roles of catalytically more 
and less active constituents in the absence of Mg are reversed in its pres- 
ence. Finally, a damped periodicity is found which indicates the pres- 
ence of two factors of an unknown sort which influence and are influ- 
enced by the inactivation of the enzyme. 


The heat inactivation of enzymes in the absence of substrate is a 
well known phenomenon. So also is the decline in the rate of enzyme 
reactions during the course of this reaction. Shutz’s law takes ac- 
count of this decline but is purely empirical and postulates no mech- 
anism. Attempts to derive this law, or approximations to it have in- 
volved diverse assumptions, including the notion of an enzyme-inhi- 
bitor complex (Arrhenius, 1907) and slow diffusion of the substrate, 
(Moelwyn-Hughes, 1936). 

In some sense, whether by mechanical blocking (Moelwyn- 
Hughes), or whatever, the enzyme seems to be undergoing an inacti- 
vation during the course of the reaction, and this process is distin- 
guished—and presumably different—from the process of heat inacti- 
vation undergone in the absence of the substrate. It is reasonable to 
suppose that a comparison of the two processes for a given enzyme 
should reveal some properties of this enzyme not apparent from a 


- consideration of one of these processes. alone. 


Hydrolysis of $-glycerophosphate by a group of phosphates has 
been studied for the purpose suggested. In general, the substrate con- 
centration has been made sufficiently high to justify our neglecting 
this as a significant variable. Also, the concentration of the enzyme 
extract has been taken low enough so that the reaction rate is simply 


83 


84 KINETICS OF ENZYME INACTIVATION 


proportional to this concentration. In all cases the solution has been 
strictly buffered at a pH of 9.1, with a maximum variation of .02, and 
the temperature has been 37° with maximum variation of .2. The pro- 
cedure has been to incubate the extract in the absence of the substrate, 
for a period of time, and to measure in gammas the amount of inor- 
ganic phosphates produced in a standard solution at varying periods 
after combining with the substrate. Throughout, we shall denote by 
p,(t) the amount of inorganic phosphates produced after ¢ hours of 
the reaction catalyzed by an extract which had been preincubated for 
7 hours. A detailed report of the experiments will be published else- 
where. 

A priori the extract might consist of any number of constituent 
enzymes, X, Y,---, acting independently. Actually, to interpret the 
results obtained, the simplest supposition seemed to be that all the 
extracts studied contained exactly two such constituents, with marked- 
ly different behavior and in markedly different proportions. In only 
one case was the presence of a second constituent doubtful, though in 
no case did it reveal itself with any certainty until the reaction had 
run for a period of 24 hours or longer. 

If there were one enzyme alone, on the proportionality assump- 
tion we could write 

dp 

—=y 

dt 
where x represents the product of the number of moles of enzyme 
present by the number of moles of substrate decomposed in unit time 
by one mole of X. Since we have no means of separating these two 
factors, we shall speak loosely of « as measuring the “amount” of X. 


If two enzyme constituents are acting independently, we have in place 
of the above equation, 


—=axty. (1) 


If X and Y are inactivated monomolecularly by heat at rates h and 
k , we have 
dx/d+ = SiMe dy/dr = ifn (2) 
whence ‘ 
Lr=%oe™, Yr=Yew, (3) 


If the rates of inactivation during the reaction are a and B, the 
process being again monomolecular, then 


dx/dt=—axr, dy/dt=-— fy, (4) 


ALSTON S. HOUSEHOLDER AND GEORGE GOMORI 85 


and, in view of equation (1), 


p7(t) =a,-(1 —e%) /a + y-(1 — €8') /8. (5) 
In case f is small, the last equation becomes 
0, (t) st 4G) a CeO Ustie (6) 


These equations result from the simplest suppositions we can make, 
and we proceed to consider the results in the case of one extract for 
which we have the most abundant data. Some of the observations 
were repeated after the enzyme had been stored in the refrigerator 
for a period of some months, and the behavior was seen to change 
somewhat during this time. 

Data from the fresh enzyme, with 7 = 0 and 7 — 24, were found 
to follow equation (6) with 


a— 0.0 nour 2 bh — & — 0.029 hour; -4,/o,— 160; 7, = 3.3. 


Here and throughout, « , and k are measured in inverse hours, while 
», x and y are in arbitrary units. The dimension of » is that of a 
mass, the dimension of x and y—mass per unit volume. Measurements 
of » were made in this as in all other cases, for t = 1, 2, 4, 8, 12, 24, 
48, 72 and 96 hours. Since 


0.029 = (log.2) /24 , 


these values of h and k amount, according to equation (3), to setting 
Gee of 2, Yor —Yo/ 2, SO that, by-equation (6): 


Do(t) =2 Dos (Ct). 


: 4 8 12 24 48 72 96t 
FIGURE 1 
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Thus the heat inactivation of the X and Y constituents proceeded at 
the same rate, whereas during the reaction, X broke down at about 
30 times this rate, Y not to any measurable degree. 

After being stored for a time in the refrigerator, these experi- 
ments were repeated and others performed, incubation being allowed 
to run also for 48, 72, 96 and 120 hours (Figure 1). All curves satis- 
fied equation (6), and except for p, , it was found that 


a=0.8hour?; h=0.00723 hour?; k=.00349 hour’; 
tof 100 3) Ya 2.2 
However, the parameters required for the ».-curve—which was run 
with half the amount of enzyme extract—were 
oS hour“ 6/0: = 96 molbem er, —— oe 
It appears, therefore, that with the extract thus stored: 


(a) During incubation a transient process takes place which is over 
fairly soon, certainly within 24 hours. 


(bo) After the termination of the transient process, the inactivations 
of both X and Y proceed monomolecularly, but at a much lower rate, 
(about 4 and 4 respectively for X and Y), than was the case for the 
fresh preparation. 


(c) The rate of inactivation of X in the presence of the substrate 
is but slightly affected. 


On the basis of available data, nothing can be said as to the nature or 
course of the “transient process,” although we may note that for p, 


———* WITH Mg 
—— WITHOUT _Mg 
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the value of x) is about twice, that of y, nearly half again the values 
to be expected on the basis of other curves. 

After a still longer period of storage, measurements of p;(1) 
alone were made for 7 = 0, 1, 2, 4, 8, 24, 48, 72, and 96 (Figure 2). 
From the graph and from equation (1) it follows approximately 


pr(1) =dp-(0) /dt=4x-(0) + y7(0). 
It was found that for 7 > 0, 
#—.0d14 hour = k= 08 2, 111.8*. 4, = 3.2. 


The value of (1), however, instead of being 115 as would be ex- 
pected from the other data, was only 94. Thus again the occurrence 
of a transient process during incubation is indicated, and here it is 
clear that this is at an end within one hour. Note that after the long 
period of storage, Y seems to have become completely protected from 
heat inactivation, whereas X breaks down even more rapidly than in 
the fresh preparation. From these measurements alone, a could not 
be determined, but the entire curve p, was run, from which it ap- 
peared that a had fallen to the value 0.12. Thus after the long stor- 
age, whereas X has become much more subject to inactivation by heat 
alone, it has become more resistant to inactivation during the reac- 
tion. 

A similar set of measurements of p;(1) was made at the same 
time as that just discussed, but with Mg added (q Mg = 2.3) with 
the substrate (Figure 2). For this it was found that 


h= 0.0525 hour; k=0; #=—151.2; y¥=3.8. 


The difference between the two values of h cannot be regarded as sig- 
nificant, and one would suppose them to be the same. The effect of the 
Mg here seems to be to render Y slightly, X considerably more active. 

Data on other extracts are considerably less complete, but serve 
to exhibit a definitely new feature, in that a itself is altered by incu- 
bation. Thus in equation (6), we must write everywhere a, for a. 
For one extract with incubation periods of 0, 24, and 96 hours, it was 
found that 


d= 0.22hour+; o.4=1.18hour*; a. = 0.69 hour -;: 
ee 00 E- 2 rt atz/ Or = 66 e--°1787, 


Now, from equation (6), z/a is the amount of p produced in an in- 
finite time by the activity of X alone, and it is this rather than x 
which undergoes the monomolecular decrease during incubation. It 
thus appears that each X molecule is limited as to the number of 
molecules of substrate that it is capable of decomposing, and that it 
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is this “capacity” factor which is affected by the incubation. The a, 
on this interpretation, measures only indirectly the rate of break- 
down of X in the presence of the substrate. Directly it measures the 
rate at which the X decomposes the substrate, and after the X mole- 
cule has decomposed its quota of substrate molecules, whether this be 
soon or late, the X molecule is capable of no further catalytic action. 

This rate factor, a, is seen to vary independently of the capacity 
factor, first increasing and then decreasing, in a manner which can- 
not be further specified in the light of available data. 

Returning, for the moment, to the extract previously discussed, 
it may be noted that since the a there was unaffected by incubation, 
we have no assurance that it was not actually the capacity factor of 
this X also which was diminished by incubation. However, in both 
extracts, it is clear that the Y is not so limited to any noticeable ex- 
tent. 

Returning now to the second extract, it was found that when Mg 
was added, the same loss in the capacity factor for the X was exhi- 
bited again, just as without Mg: 


27/ dr = 808 €--02787, 


The a’s, and here also the y’s, showed the same general type of varia- 
tion as the a’s without the Mg: 


Yo = 0.64 aes en ee een ta 
Q) = 0.0556 hour?; og,—0.087 hour?; a ,—0.07 hour . 


On a third extract, only four curves are available, for 0 and 24 


hour incubation, the reaction run with and without Mg, but they show 
striking features. We write the equations completely: 


Po= 1.76¢ + 114(1 — €%-07#) | 

Dos = 1.26¢ + 80(1 — e-24) , 
without Mg, and 

Po = 1.58t + 213 (1 — e068) | 

Do, = 1.58¢ + 150 (1 — e-°-047#), 
with Mg. 

The values here obtained are extremely interesting. If we de- 
note by a bar (e.g. x) the values obtained with Mg, we note that the 
logarithms of the ratios of the y’s, the (%/a)’s, the (x/a)’s, and the 
a’s are all very nearly equal to 0.014 or 0. 015 times 24, indicating 
monomolecular breakdown with time-constant about 0.015 inverse 
hours. Hence the same is true of the x. On the other hand 


Ly = Ag (Xo/O%) =7.98 3 t= 8; 
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Yo = Yos = 1.58. 
Thus the heat inactivation of the Y seems to be completely reversed 
by Mg. As for X , it seems to be subject to several processes during 
incubation. 
Neglecting, for the moment, the effect of Y, which is, in any 
case, purely additive, let us write 


dé/dt =— 1 dp/dt, 
dp/dt =n, we 


where & represents the number of moles of X , u the number of moles 
of substrate decomposed in unit time by one mole of X, and 4 is the 
number of molecules of substrate that each molecule of X is capable 
of decomposing. Then we find 


P—so Atl — er), (8) 


Here &, is the number of moles of X present at the beginning of the 
reaction. If we consider the non-magnesium curves, we find that dur- 
ing incubation, apparently, & 4 decreases, u// increases, & uw remains 
constant. Evidently, therefore, 4 decreases, & and u remaining un- 
affected by the incubation. But from the magnesium curves, appar- 
ently, w/A_and & both decrease in the same way whence, presumably, 
4 remains constant, & and uw both decrease. 

It appears, therefore, that the modes of catalytic action of both 
constituents X and Y exhibited in the presence and absence of Mg 
are distinctly different. Thus & and wu both show declines with, but 
not without Mg, while the “capacity” factor, 2, shows a decline 
without but not with Mg. Whatever the change undergone by X dur- 
ing incubation, it amounts to a limitation only of the capacity factor 
A without Mg, but with Mg it diminishes both the amount available 
for catalysis, and the rate at which it acts. 


We have been designating, in general, by X the most active, by 
Y the least active enzyme constituent, though without any intention 
of implying identity of the various X’s and the various Y’s. In the 
present instance we have not even any assurance that the X of the 
no-magnesium reaction might not be identical with the Y of the mag- 
nesium reaction, and vice versa. There is, in fact, a suggestion that 
they may be, since, as noted above, in the one case it is the x, in the 
other case the y which remains constant. On this assumption the X , 
say, which does not break down during incubation, acts rapidly but 
with heat-labile capacity factor in the absence of Mg, while in the 
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presence of Mg it acts only very slowly and with a capacity factor 
which could not, in any case, be detected (u being too small). On the 
other hand Y, which decomposes monomolecularly during incubation, 
acts rapidly in the presence of Mg, but at a rate w which also de- 
creases during incubation, whereas in the absence of Mg, its wu is too 
small to be detected in the exponential. 

Nothing more definite than this can be said, and further specu- 
lation would hardly be profitable here. It is intended only to outline 
such conclusions as can be drawn from the available data, and to sug- 
gest a few possible developments. We mention, finally, one other ex- 
tract upon which further work is planned. When the graph of p;(1) 
(with Mg) against 7 is drawn on semi-logarithmic paper, one recog- 
nizes a curve of the following type: 


p,(1) = ae — be-fbt sin yr. (9) 


Such a curve can be obtained only if the system contains at least three 
constituents, and the simplest form of the equations is 


dt/dn==— Ad — Bw — Cz, 
dw/dr=x, (10) 
dz/dz=w 
Hence, assuming only a single actual catalytic element X , this X be- 
sides catalyzing the decomposition of the substrate, also “catalyzes” 
the development of a substance or state or condition W , which, in 
turn, contributes to the development of a substance, state, or condi- 
tion Z. On the other hand, both W and Z act to break down the X. 
Nothing can be said as to the physical character of W and Z without 
physical or chemical measurements of a different sort. They could 
be any deviations of physical or chemical parameters from some “nor- 
mal” equilibrium values, and as such, w and z could be either positive 
or negative. 
The authors are indebted to Dr. H. D. Landahl for numerous 
helpful suggestions. 
This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 


LITERATURE 


Arrhenius, S. 1907. Immunochemistry. Leipzig. 


Moelwyn-Hughes, E. A. 1937. “The Kinetics of Enzyme Reactions.” Erg. d. 
Enzymforschung, 6, 23-46. 


BULLETIN OF 
MATHEMATICAL BIOPHYSICS 
VOLUME 5, 1943 


NOTE ON ANODAL EXCITATION 
ALSTON S. HOUSEHOLDER 


THE UNIVERSITY OF CHICAGO 


A general, necessary and sufficient condition for the occurrence of 
anodal excitation is given on the basis of Rashevsky’s two-factor theory 
of excitation, expressions for the anodal rheobase and the minimal shock 
intensity following as simple special cases. 


Rashevsky’s equations (Rashevsky, 1938) are equivalent to 


—=KI— ke, 
dt 
ai (1) 
== MPa, 
dt 
with excitation occurring when 
o=e-—-72l, (2) 


since it is no restriction to so choose the units that the threshold 
(= 7. — & in Rashevsky’s notation) is equal to unity (cf. Young, 
1941). 

For brevity, we shall speak of “anodal” excitation in general as 
that occurring only after the cessation of the stimulus. It is a simple 
matter to show that if at any time when o < 1, and either « > 0, or 
j > 0, the stimulus is discontinued (cathodal break), excitation can- 
not ensue without further stimulation when k > m (Young, 1941), 
but that the stimulus being discontinued wheno <1,2e<0,7< 0 
(anodal break), excitation may ensue. If, at break or any time there- 
after, the « and 7 have values s, < 0, and 7, < 0, the “excitatory 
state’ o decays according to 


car e-kt — do emt, (3) 


and for excitation to follow t seconds thereafter, it is necessary that 


at this time 
a= & Code aed0 emi — 1 ’ 
s—mj,e" —ke, e220, 


or, what amounts to the same, that at this or some later time ft, 
OL 


92 ANODAL EXCITATION 


Gale o =—0. (4) 


Consider the relations between ¢) and 7, implied by these latter condi- 
tions. If we solve the equation (4) for ¢, and substitute into the in- 
equality (4), we obtain, after elementary transformations, 


(—Jo)* ae kk 


(Ct sy m”" (k = mn) 2 


(5) 


as the necessary and sufficient condition for anodal excitation. 
In Rashevsky’s theory, k 2 m. To examine this condition more 


closely set 
k= (1 +a)m, OsaL0e (6) 


The fraction on the right then becomes 
[(1 + a.) 1+a ppd Veer 


and the relation may be expressed, after taking the m-th root of both 
sides, 
(Jo) te 6g CLE a) en (7) 


lrewerurthernctennat 
(i= eo) 8? o* = (To) 2) ee 


both factors approaching unity in the limit as a — 0, we may take the 
logarithm of both sides and write 

(1 + a) log(—Jjo) 2 log (—eo) 
+ [(1+ a) log(1 +a) —aloga], a) 


where the quantity in brackets is positive and, in fact, increases with 
Qa. 
For convenience, let 


= 102 (2 Jo) a — 10n aea (9) 
If we take the equality sign, we have a family of straight lines 
(l+a)y=e2+[(1 +a) log(1+a) —alogal], (10) 
and these lines envelope (all are tangent to) the curve 
e—e=1, (11) 


the point of tangency of the line defined by any « being given by 


k 
tae Ors ee, : Ui NOS bilan loa eA 
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The curve (11) has for its asymptotes the x-axis and the line y = a; 
it rises monotonically for all x , and has the y-intercept log 2. 


The curve (11) is, of course, merely the logarithmic graph of 
the line 


Melo) an Geeo)e—— lv. (13) 
Hence, in terms of the ordinary graphs, we may say that the curve 
eos (ey) (Loe) (14) 
is tangent to the line (13) at the point given by 
—eQ=m/(k-—m) , —j=k/(k—m) (15) 


and the condition for anodal excitation is that the initial point (—«), 
-(—jo) shall be in the upper region bounded by the curve (14) and 
the line (138). In terms of the logarithmic graph, the region is defined 
by the line (10) and the curve (11). 

The expression for the anodal rheobase is easily derived from 
(14) if we set therein 


—e=KR/k , -JFo=MR/m, 


the asymptotic values of « and 7 produced by a steady current —R,. 
We obtain the form 
k m Km (km) 
R= ; (16) 
Hes Oye ME) 


somewhat different in appearance from that given by Rashevsky 
(1938), but easily seen to reduce to the same. In like manner, if we 
set 

—& = K Qu ’ —jo=MQz., 


where Q, represents the minimal instantaneous quantity of current 
capable of producing anodal excitation, we obtain 


Km km) gki(km) 


Qa = (17) 


= M/k) mn (k-m) k —m : 


In particular, since the minimal cathodal shock required for excita- 
tion is (K — M)-— (Young, 1941), this is less than Q, if, and only if, 
(l-- A)AS> a? (lo) oO (18) 
where 
earl eA ae aA 0G (19) 


We may note, finally, that by combining a steady current with 
an immediately succeeding shock, ¢ and j can, in principle, be varied 
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independently of one another, since with such a combination 


e=K(1+kQ)/k, 
j=M(I+ mQ)/m 


and, m being different from k, any assigned « and 7 can be produced 
by suitable choice of J and Q. These relations, together with equation 
(5) or (7), make possible the testing of the theory by means of anodal 
stimulation over a somewhat wider range than is provided by (16) 
and (17) alone. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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ON THE THEORY OF AMOEBOID MOVEMENTS 


N. RASHEVSKY 
THE UNIVERSITY OF CHICAGO 


It is shown how the application of the principle of maximum energy 
transformation, discussed in previous papers, leads to expressions which 
give the probability of occurrence of a pseudopod of a given size and 
duration in function of other physical constants of the amoeba. 


In previous papers (Rashevsky, 1939; Stanton, 1943), it has 
been suggested that the form changes of an amoeba may be approxi- 
mately treated as a quasiperiodic process of elongation and contrac- 
tion. The actual irregularities due to simultaneous protrusion of sev- 
eral pseudopods in different directions would be considered as due to 
disturbing factors to be accounted for by a second approximation. 

In treating quasiperiodic changes of cellular form from the point 
of view of the dynamics of the diffusion forces, one might think of 
applying the equations of deformation to each individual pseudopod 
in the same way as they are applied to the cell as a whole. It is, how- 
ever, necessary to introduce the restraint of a constant total volume, 
and this makes the whole problem mathematically exceedingly diffi- 
cult. 

If, however, we treat the quasiperiodic changes from the point 
of view of the principle of maximum energy exchange (Rashevsky, 
1943a, b), then the solution of the problem becomes in principle rath- 
er simple. 

Let 

N(v,1,0,¢) dv dl dé d¢ (1) 


be the relative frequency of formation of pseudopod whose volume 
v, length 1, duration 6 and direction ¢, as given by the angle in a 
plane from a fixed direction, lie correspondingly in the intervals v , 
y+dv;l,l+dl;6,6+d0;and¢,¢+d¢. 

If N, is the number of pseudopods formed during a sufficiently 
long time, then expression (1) represents the ratio of the number of 
pseudopods whose v , 1, 9, and ¢ lies in the specified intervals to the 
total number N,. Hence 


$2 Se Se S27 N(v,l, 6, 4) dv dl do dp=1. (2) 


The protrusion of a pseudopod influences the ingestion of food, 
95 
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and the ingestion of food constitutes the positive energy flow into the 
amoeba. Depending on the different assumptions which we may make 
about the mechanism of food ingestion, we shall arrive at different 
expressions for the amount of food ingested by a pseudopod char- 
acterized by given values of v,1,6, and ¢. In general, we could ex- 
pect the amount of food ingested to increase with v, l, and 6. For 
instance, we may assume it to be proportional to the surface of the 
pseudopod. Since the crosssection area of the latter is v/l, therefore 
the surface will vary like /\/v/I and the amount of food ingested will 
vary like 6l\/v/l. Similarly a number of other assumptions may be 
made. 

If, in general, the concentration of food is nonuniform, the ex- 
pression for the amount of ingested food will involve ¢. 

For every particular assumption which we make about the pro- 
cess of ingestion, we can in principle determine the amount f,( v, 1, 
6, ¢) of energy flowing into the amoeba during the existence of a 
given pseudopod. 

But the protrusion of a pseudopod, by increasing the total sur- 
face, will increase the oxygen concentration in the cell. If we assume 
that the respiratory reactions in an amoeba are of the kind studied 
by H. D. Landahl] (19389), then an increase in oxygen consumption 
will result in an increase in the oxidation rates, and hence in an in- 
creased loss of energy, due to the exothermic reactions, in which the 
energy is dissipated in the form of heat. By the usual approximation 
method, the change in oxygen concentration for a given pseudopod 
may be calculated. A procedure similar to that used by H. E. Stanton 
(1942) may be employed. 

In this way we may calculate the amount f.( v,1, 6) of the nega- 
tive energy flow caused by a pseudopod. The net positive energy flow 
into the cell is f= f, — fe. 

We shall postulate that the distribution function (1) is such as 
to maximize the average total positive energy flow, due to all kinds 
of pseudopods. 

Let us consider for simplicity, first, the case when v, l, 6, d 
assume a discreet set of values 0; , lk, Om,¢r,4,k,m,7T=1,2,---n. 
Let the frequency of occurrence of v;, I: , 9m, ¢r, be N 

The total energy exchange is given by 


ikmr « 


a Niame FU; lis Om 5 br) . (3) 


i,k,m,r 


If we denote by v ,1, @ and ¢ the constant average values of the vol- 
ume, length, duration and direction angle of a pseudopod, then ex- 
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pression (3) is to be maximized subject to the conditions 


SNe v; DEN thane Gk i 


ikmr ikmr 
>> N ike On = 6; > Nixme br =e $. oe 
ikmr ikmr 


Denoting by 1, u, », and w Lagrangian multipliers, the above re- 
quirement gives 


f(v; A L, Om 5 hr) ar AV; oe ub, Ar VOm wor =). (5) 


The Nixm,’s drop out entirely. Hence, with those assumptions the 
frequencies of occurrence of the different pseudopods are not intrin- 
sically determined. They may be assumed to be due to pure chance, 
and N(v ,1, 6) would be a normal distribution, independent of ¢. 

Equation (5), which for the continuous case becomes 


f(v,t,6,¢6) +aAv+ pl +76 + wod=—0, (6) 


can in principle be verified experimentally, since it holds for each 
pseudopod, and the parameters 1, u, » and w can be determined from 
observation. 

Though the distribution may not depend explicitly on ¢, yet that 
quantity affects the situation through relation (6). While a pseudo- 
pod in any direction may be equally probable a priori, yet for the 
pseudopods in different directions, the relations between v, 1 and 6 
will be different, for they will correspond to different values of ¢, 
which now plays the role of a parameter. As remarked before, the 
dependence of ¢ will enter only through the nonhomogeneity of the 
external medium. The equal probability for pseudopods in any direc- 
tion would result in a zero average displacement of the amoeba as a 
whole. This might be the case for an absolutely homogeneous external 
medium, for the amoeba moves towards regions of higher food con- 
centrations. 

We may, however, consider a more general case in which the dis- 
tribution function is itself intrinsically determined by the require- 
ments of a maximum energy interchange. 

When a pseudopod is protruded, this increases the positive energy 
exchange because of increased food ingestion. Remembering that an 
amoeba ingests rather large particles, the distribution of which in the 
medium is not necessarily uniform, we see that if a pseudopod of the 
same size, shape and direction is protruded too frequently, the aver- 
age increase of food ingestion per protrusion will decrease because of 
exhaustion of food in that region. Hence the function f, will depend 
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not only on v, 1, 6 and ¢, but also on N(v, l, 6, ¢). Considering 
again the case of discrete values, we have 


Fi = Fi. Urs Oo ee aoe (7) 
Instead of expression (3) we now must maximize the expression 
>, IN itens [7 (%, kee Ams or > Nikmr) ay fo(vi, Lk. , Om) J ’ (8) 
ikmr 
subject to conditions (4). We now find 
of 
ipa Gee ke; Gia dr» Nixmr) =~ fe(vi, L, 5 Om) ay Ninmr ae 


= (9) 

ae VDE SF ub, SEM Sar Od, =0. 

This is an algebraic equation in Nizm,, which determines the lat- 

ter in function of vi, k, 0m, ¢r,4, uw, » and w; thus giving for the 

continuous case the function N(v ,1,0,¢,4,u,¥,@). The paramet- 

ers 1, u,v and w are determined from equations (4), which for the 
continuous case become: 


fs fON(0,1,0,6,4, 0,0, dv dl do dp = (10) 


and three corresponding equations for], 6 and 6. We have four equa- 
tions for the determination of the four quantities 1, uw,» anda. 
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The new fundamental equation, based on the principle of maximum 
energy exchange, is applied to the case of elongation and constriction 
of a cell. A simple case is treated as an illustration, and it is shown 
that the elongation curves thus obtained are of the same character as 
in the old theory. 


Presupposing familiarity of the reader with previous publica- 
tions (Rashevsky, 1943 a, b, c, d; hereinafter referred to as [a], [b], 
[c], [d]), we shall apply the fundamental equation of [c]: 


dx; oF oH 


a 1 
J dt ae Ox; a) 


to the problem of elongation and constriction of dividing cell. 

Neglecting the energy dissipation in the diffusion field itself [d], 
the expression for F in a cell is some function of the concentration ¢; 
of the different metabolites. Those ¢, are in their turn determined 
by the different diffusion coefficients, permeabilities, reaction rates, 
as well as by the two coordinates, which roughly describe the shape 
of the cell, namely its length 27, and the width of the neck 7, ([b], 
[d]). The function F'(7,, 7) is determined in any particular case by 
the approximate solution of the diffusion equation. 

Considerations of viscous flow determine the quantities f,, and f, 
which give us the viscous forces f,,(d7,/dt) and f,(dr/dt) of resist- 
ance to the elongation and constriction. H is given by the potential 
energy of the surface tension. If S denotes the surface of the cell, 
then 

H=yS. (2) 


The question may be raised as to whether the potential energy 
due to the diffusion forces should also be included into H. Although 
the diffusion forces are formally derivable from a potential, yet their 
very nature is based on non-conservative phenomena, and therefore 
such energy considerations as have been used before (Rashevsky, 
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1938, Chapter viii) have always been felt as somewhat uncertain. 
There is also this fundamental limitation of the diffusion forces: their 
effects vanish with infinite D,. This is extremely implausible bio- 
logically, and has been always the weak point of the theory. In cases 
of protoplasmic streamings accompanying division, there remains the 
old difficulty ; namely, although the streaming itself is rather an indi- 
cation of the existence of the forces, their distribution may be upset 
considerably by the streaming, so as to strongly modify their effects. 
It is more likely therefore that the diffusion forces should not be gen- 
erally included into H, although they might be included in special 
cases. We shall therefore omit them here. Other sources of potential 
energy, such as possible elastic forces should, however, when present, 
be included in H. 

The system of differential equations which describe now the 
elongation and constriction is: 


dr, s or (%4 ’ r) oS (r1 ’ r) 


ita puaes Leet ar, 
dr OF (71,7) oS (11,7) (3) 
iiamaebigay haem me 


As an illustration merely, we shall apply equations (3) to the 
simple idealized case when a cell elongates without constriction, so 
that the second equation drops out, and when we have only one sub- 
stance, produced exothermically at a rate proportional to its concen- 
tration. We shall also for simplicity put D. = h = ow, and denote D; 
by D. 

The frictional force may be computed by considering the elongat- 
ing cell as a plastic rod of length 7, and width r.. For such a rod we 
have 


Ld? sh Le 
r, dt 3n8” ee 
where P is the total force, and s the crossection. 
Remembering that 
— To ; 3 Tatas — MEL = Vaz T1727 == Te ; (5) 


solving equation (4) for P and ‘lifiinatite Yr, by means of (5 
find for the force the expression } (5), we 
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dr, ae 9nV dr, 


14 os 3 6 
hee VCE yh a (6) 


Denoting by w the heat of reaction per gram of substance pro- 
duced, and remembering that the reaction is exothermic, we have, 
denoting the rate of reaction by ac: 


F=—-Vwac. (7) 
But 
Co Ty? 12? 


LESS SoS era (8) 


l—an = 8D (2r2 +72) 


— 


Hence, eliminating r, by means of expression (5), we have 


27,” +7o3/r 
F=- Ar ae wt le (9) 


ee 
BD(2r lis ae aR, 
1 
and therefore 


dF 9V* Dwa? ¢ (I Mae Or dy 


eee Whether oe ee 


L 3 
[3D (2r,2 + =) — aren]? 


For r, > 7, dF /dr, > 0. It is zero for 7; = 7 and for 7, = o. 
Introducing (6), (10) and (2) into the first equation (3) we 
have after rearrangements: 


dr, 2V Dwa? Cot bie a sisal Was AT Y dS 

dt ne, Li) To? 

Lol (2¢2 -- ae ine 71 |< 
1 


The last term in equation (11) may be calculated by approximating 
the cell by an ellipsoid. The general expression is rather cumbersome. 
For small values of elongation 7, — 7, we have (Landahl, 1942) 


Y 
S=4ar? + 1-6 ard (= — 1)? (12) 
and therefore 
aS 
—=1-6a(1-—71). (13) 
dr, 


Since for very small values of 7, — 7 the first term of the right 
side of equation (11) may be written 
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6VDwa? Co To° 7 
——_————_ (711-1), 
my (9Dr,? a, Oro) 2 


the requirement that dr,/dt > 0 for 7, = 7 leads to, remembering 
equation (138), 
8D we? C) To° t gt 
(ODT — G734)- 3 To 


(14) 


which is satisfied only when 7, exceeds a critical radius, determined 
by making the two sides of expression (14) equal. The universal con- 
stant z is one of the factors determining the critical radius. 

For very large values of 7, — 7) considering the cell as a very 
long cylinder, the surfaces of whose ends can be neglected, we have 
S = 2ar,’ \/r, and therefore, dS/dr, = 2ar,°/?/7,". The last term of 
equation (11) now becomes (2y/3”) (7,°/2/7,*/7). Since the first term 
of the right side of equation (11) tends to a constant value for large 
values of 7, — 7, therefore when a certain value of 7, is exceeded, 
dr,/dt becomes zero, even though inequality (14) is satisfied. Hence 
the general shape of the elongation curve is the same as in the old 
theory. Elongation however occurs even for D, = h= o. It occurs 
also for D; = D., = wo, h—finite, when no volume diffusion forces can 
be present. . 

From the comparison of calculated and observed elongation 
curves we may determine the constant r. It can also be determined in 
principle from the critical radii, although this would not be practical. 

The author is indebted to Dr. Herbert D. Landahl for verifying 
the calculations. 
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The theory dealing with recognition learning is generalized to in- 
clude five variables. These are the errors, trials, possible number of 
choices, the number of items and the fraction of prompting. The theory 
is then extended to include recall learning. This extension necessitates 
the introduction of a variable which is the number of trials without 
responses, i.e., the third category of judgment. The strengths of reward 
and punishment are also introduced as variables. Thus a curve in nine 
dimensions is determined. 


In previous studies we have considered a neural mechanism which 
exhibits properties of discrimination and conditioning (Landahl, 
1941a, b, hereinafter referred to as I and II respectively). It is our 
purpose here to generalize the previous results for the case of recog- 
nition learning as well as to introduce recall learning. We presuppose 
familiarity of the reader with the previous two papers. 

1. Recognition learning. We consider first the situation in which 
recognition of the learned response is all that is required. The experi- 
mental situation to be considered is much the same as that described 
in II. Let there be M pairs of stimuli whose association is to be 
learned. Let one of each pair, S;(4 = 1,--- M) be the key stimulus 
which when presented will, after sufficient learning, illicit the re- 
sponse normally given to the other stimulus of the pair. This latter, 
S’';(g =1,--- 'N), is one of a group of N stimuli arbitrarily chosen, 
except that the group includes the stimulus to be associated with the 
key stimulus. When a stimulus S; is presented, the subject is to choose 
one of the S’; stimuli. Thus he has N possible choices. We assume that 
enough time is allowed so that a choice is made on each trial n. If 
the choice is correct, the subject is so informed or rewarded, and this 
contributes to c;, the cumulative correct responses when S; is pre- 
sented. As we consider the M stimuli S; to be of the same intensity 
we may drop the subscript. Thus ¢ is to be identified with the average 
of the experimental c values. If the choice is wrong, it contributes to 
the cumulative errors, 2. The subject is then informed of the error 
or punished. He may or may not be shown the correct response sub- 
sequently. In general, the experimenter may tell the subject as to 
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which is the correct response a certain fraction (1 — f) of the times 
in a particular experiment. If f = 1, there is no prompting, and if 
f =0, prompting follows each trial. Thus the variable f must be on 
the closed interval from 0 tol. 

We consider that the neural mechanism has essentially the prop- 
erties discussed in I and II. We assume that it is sufficiently complex 
so that conditioning may take place between any S; and S’;. The con- 
ditioning may be such that the tendency for S; to produce the response 
to S’; is either strengthened or weakened. For simplicity we shall neg- 
lect any bias so that at the first trial the N stimuli, S’; , are equivalent. 
Let «) be the initial value of the excitation factor corresponding to 
any one of the N stimuli. After c correct responses to the S’; stimulus 
associated with one particular S; stimulus, the effective value of the 
excitation factor for the correct stimulus will have increased to some 
value e, due to the conditioning. The N — 1 wrong responses will each 
have decreased by various amounts. But, as in II, we shall not con- 
sider the individual wrong responses, only their total. Thus the aver- 
age value of the excitation factor for the N — 1 wrong stimuli is ail 
we shall use, and it is given by some value e, after w wrong responses 
have been made. 

As only one of the S’; stimuli can be chosen, we again consider 
that there are cross inhibitory fibers which tend to block out all but 
the stimulus which produces the strongest excitation. Introducing 
the threshold fluctuations, we shall expect each stimulus S’; to have 
a certain probability P; of being the stimulus to pass through, and 
this probability will be greater for the stronger excitations. In gen- 
eral, this relation will depend upon each of the stimuli, upon the con- 
ditioning which has taken piace in all possible combinations, and 
upon the number N . We shall, however, use the approximate expres- 
sion (I, eq. 9) 


w 


lo +k(e.— &») =0, 1 
Bay ( (1) 


where k is a constant and P,, is the probability of a wrong response 
at m trials. There being only two categories of judgment, the prob- 
ability of a correct response, P,, is equal to 1 — P,,. 

Considering again that the strength of conditioning varies lin- 
early with the number of correct responses, we will have for «, an 
amount be + b(1 — f)w added to e), where b is a constant. In gen- 
eral we might have introduced two constants, as the guessed response 
and the response due to prompting are not identical. The mechanism 
does not allow this distinction without further complications. We 
shall thus consider only the simpler case. However, it may be pointed 
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out that by permitting f to have values somewhat different from 
their preassigned values, this effect can be included. Thus, if f must 
be assigned a value f’ instead of f , we would interpret the discrepancy 
to mean that the ratio of the actual coefficient of (1 — f)w in (2) is 
given by (1 — f)/(i — f’). This ratio could then be introduced as an 
additional parameter. Subsequently we shall introduce factors which 
will complicate such an interpretation, so that the above would only 
be an approximation. 

The act of responding to any stimulus, together with the subse- 
quent punishment or reward will, in general, react on the state of the 
conditioning at any center. In terms of the neural mechanism, this 
means that inhibitory fibers leading to various conditioning centers 
may be excited, thus tending to disrupt the conditioning. After n 
trials, there will have been Mc correct responses altogether. Thus an 
amount roughly proportional to (M — 1)c is subtracted from the con- 
ditioning at a center corresponding to a particular stimulus S; . With- 
out loss of generality, we let this amount be proportional to Mc. In 
general, one would expect that the inhibitory effect would increase 
with Mc at a decreasing rate. Making such an assumption increases 
the complexity considerably. However, if instead of M, we write a 
function of M which increases at a decreasing rate, no great difficulty 
is introduced. Thus we may write for the amount subtracted; ¢(M)c, 
which to a first approximation may be written ¢’Mc, or to a better de- 
gree of approximation ¢’(1 — e?")c/u, ¢ being a suitable function of 
M , and ¢' and wu being constants. In a similar manner we must sub- 
tract a quantity 7(M)w for the effect of the wrong responses, 7 (M) 
being similar to ¢(M). Furthermore, each S’; stimulus is itself re- 
sponded to in some way, although not recorded by the observer. Thus 
one should subtract a quantity 6(MN)7 as this type of response takes 
place MN times at each trial. 6(MN) is a function similar to ¢(M), 
but probably very much smaller in value. On adding together the 
various terms algebraically , we may obtain an expression for ¢é¢. 

An expression for ¢, can be obtained analogously. The princi- 
pal term in the contribution to «, due to the conditioning to avoid 
wrong responses is 8 w/(N — 1), where £ is a constant analogous to 
b. For each wrong response we have an amount £ at one of the cen- 
ters corresponding to one of the (N—1) wrong responses. Thus an 
amount £/(N — 1) is contributed to the average for each of the w 
wrong responses. This result may also be obtained if we recall that 
learning to eliminate an error is 1/(N — 1) times as effective as 
learning a correct response. Similarly, each of the other terms one 
should add would be divided by (N — 1). Unless the effect of punish- 
ment is much greater than the effect of reward, we may then neglect 
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additional terms, leaving Ey = & — Bw/(N — 1). 


Subtracting the expression for é» from the expression for ¢é, 
we obtain 


Eo — Ep = (yn—dbw)/k, (2) 

where 
y=kb — kC(M) — ko(MN), (3) 
6=kbf — kB/(N —1) — kfC(M) + kn(M). (4) 


Substituting (2) into (1), and dw/dn for P, and eliminating ¢ 
by c=n-— w, we obtain a differential equation in wand n. Integrat- 
’ ing this equation, we have 

N 
w= - log MRED AOA OMNES (5) 

6 (N -"1)"6: 6-2") = CN 1) oar Nes 
Equation (5), with y and 6 replaced by their equivalents in (3) and 
(4), determines a surface in a five dimensional space defined by the 
variables w,n,f, M,and N. In some experimental situations, one 
could, perhaps, neglect the terms in f, 6, and 7, as well as the uw in 
¢(M), thereby determining very roughly the surface by two para- 
meters. This is essentially what was done in II for the special case in 
which f =—1,M=N. The term in 6 would be expected to be small 
so that 6’MN should be ample. The term in 7 could be similarly treated 
- unless the punishment is considerable. Even if so, one could perhaps 
introduce the same uw as in ¢(M). The result would still be only six 
parameters for the determination of a surface in five dimensions. 
Note that & multiplies each of the parameters in (3) and (4), and 
thus is not an additional parameter. 

In order that ¢ and w remain positive and real, and that d?c/dn? 
be positive for n > 0, it is necessary and sufficient that 


¢=y—6+6/NZ20. (6) 


This is the requirement necessary for any learning at all. If also 
y = 0 (but not both equality signs holding simultaneously), then 
w—> We. > 0 for n > o, i.e. learning is complete. But if y — 6 + 
6/N > 0 and y < 0, then w > o and dw/dn < 1 for n > o , that is, 
the learning is incomplete. The condition that there be an interval in 
M , the number of items to be learned, in which learning is incom- 
plete is 

: C 1B — (N-1) f6'] > (N-1) 7 (b—- 8), (7) 


provided that a positive, real M satisfies ¢ = 0 and if we set 6(MN) 
=6MN . From inequality (7) we find that in order to have any “in- 
complete” learning, 6 and N must be small enough. Furthermore, if 
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(7) is satisfied for f = 0, it may no longer be when f=1. That is, 
the interval for incomplete learning is diminished when prompting is 
reduced, and this is especially effective when N is increased. 

The value of N(dP./dn)/(N — 1) forn= ; 


N d?c/dn2/(N —1) =4, (8) 


gives the rate of change of the probability of a correct response cor- 
rected for chance as a function of the number of trials, and is a meas- 
ure of the rate of learning. This measure decreases monotonically 
with both M andN. 

The parameters b and ¢’ will increase with the strength of the 
reward. If the increase is roughly proportional, and if the increase 
in b is greater, the rate of learning will increase with the strength of 
the reward. If these relations are not proportional, the dependence 
of the rate of learning on reward is more complex. The parameters 
6 and y will depend on the strength of the punishment. If both are 
proportional, the rate of learning will increase with increased punish- 
ment for small enough N , but will decrease for large enough N. That 
is, in a sufficiently complex situation, increasing punishment may be 
detrimental to the learning. If the 6 and 7 are not proportional to 
the amount of punishment, the dependence of the learning upon the 
latter is more complex. 

2. Recall learning. Thus far we have considered only recogni- 
tion learning. When we treat the case of recall learning, we simply 
remove the ¢) values corresponding to the S’; stimuli since they are 
no longer presented. With successive trials, however, these stimuli 
are learned, that is, there develops an e’) which plays the same role.as 
€ , tending to produce the response which is normally illicited by the 
particular S’; stimulus. This « value is of the form [b, — ¢,(M)]n, 
the stimulus being presented on each trial. Clearly we cannot require 
a response each trial, and thus we must be dealing with the case of 
three categories of judgment, the intermediate category being failure 
to respond. In terms of the mechanism, this means that there is a 
threshold h which «. — ¢~ must exceed for a correct response to be 
given. This quantity h is not a constant in general (Landahl, 1939). 
We shall postulate that h decreases only because of the increase in 
the e«’, values. In general, h would decrease to some positive limit 
which might be zero. Unless this expression is linear, considerable 
difficulties arise in obtaining solutions. Thus we shall set ny 


h=h — bn — ¢,(M)n=hy — yin/k (9) 


when the right member is positive, otherwise h = h,. In order to fur- 
ther simplify the results and avoid introducing an additional para- 
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meter, we shall take h, = 0. 

For ¢- we have an expression similar to that of the previous sec- 
tion, except that f = 0, and that « and the term in N vanishes, that 
is @=0. But it should be noted that the terms in w arise only where 
errors are made. Let w, be the accumulated failures to respond. An 
error is defined as a response to one of the M stimuli S; which is not 
the correct response. If any other response is given, it is counted as 
a failure to respond. Thus we set M=N. 

For ¢, we have a term Bw/(M — 1) multiplied by some function 
of n which is small initially but approaches unity. We shall, however, 
neglect this variation and take ¢, to be equal to Bw/(M — 1). If B 
turned out to be somewhat smaller in this case as compared with that 
discussed above, this would be expected. One might introduce a fac- 
tor of the order of one half in this term with justification. 

We continue to use the approximation of the form of (1), except 
that now we have the three category case. Thus instead of (1), we 
have (cf. Landahl, 1939), 


w 


+k(é-— em +h) =0, (10) 


log 


M(1 — P.) = 
0 er ee ee (11) 
MEAL 


These expressions are defined for all values for which the term in 
parentheses multiplied by & are positive. Since e, — e, and h are al- 
ways positive, the first expression is defined for all values under con- 
sideration. The second expression, however, breaks down when 
Ec — & —.h <0, since for large enough negative values of this term, 
P, as defined will exceed unity. This difficulty can be avoided by in- 
troducing a new definition for P, in which continuity in P, and in the 
first derivative of P, with respect to (e, — e,) is required, as well as 
that P. should have a lower bound of zero. If an exponential type of 
curve is again used, these conditions lead to 


log MP. + (M —'1)k(h — e. + &¢) =0 (12) 


for h — e. + &» > 0. 


The value of ¢, — e~ is given by yn — dw, y and 6 being given by 
(3) and (4) but with 6 = f = 0. Now let m, be the solution of 
he eo + 8, == O-or 


ho — yn, + 6 w(n,) —w1n,=0, ; (13) 


yi being the coefficient of n in (9), w(n,) being determined from 
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(10). We shall solve only for P. and P,,, that is, dc/dn and dw/dn. 
For n = h./y; we have from (9) and (10) 


(Magen) ee 2! (yay) 


Py = dw fda —TJ—H 
M (py, — p)e — (M — 1) 8(e%P" — 1) 


(14) 


For n = n, we have from (12), eliminating w by the solution 
of (14) 


e7 (M-l)kho 
P.=dc/dn= 
M — == = —kho ( e(yi-y)n — M- 
‘x | (ny) CE) rare (ery) 1) eueDOeDA, (15) 
M(y: —-7) 


while for n, = n = h,/y:, we have from (11), again eliminating w 
by (14), 
(M — 1) (71 —= y) e- (xy) 0 ekho 


P,.=1 = Sa 16 
iM (yy) —(M — 1) 6 6 (e921 1) on 
For n 2 h,/y1, we have 
Pee ee (Mh) y er") 7 (My — (AE 1) iy 
(17) 


eee 1) (ys) 
a VEL) 6 (eet e-1") ] 


Throughout, dw,/dn is obtained from dw./dn = 1— P. — Py. 

Equations (14) to (17) define a surface in terms of the variables 
w,ce,nandM, w, being given by n — c — w. Three new parameters 
have been introduced, h, b’ and 2’; b and ¢ can be determined from 
an experiment corresponding to requirements discussed above. But 
a new variable w, , has been introduced, as well as a new experimental 
situation which is equivalent to another variable. 

It should be noted that in the case of recognition learning, we 
did not consider the possibility of failures to respond. It is possible 
to introduce a threshold h in that case also. However, in such experi- 
ments, ample time is generally given for some response to occur. Fur- 
thermore, the motivation to make a response each time is fairly 
strong. It would thus be necessary to set a time limit for each trial 
to make the experimental situation correspond with the case in which 
an h is introduced. But in the case of recall learning, it is clear that 
for the first few trials, there will be no wrong or correct responses. 
Thus we see that the relationship between recognition and recall 
learning has something in common with the relationship between the 
two and three category cases in discrimination. 
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Although not explicitly introduced thus far, we have two addi- 
tional variables. These are the strength of the reward, 7, and the 
strength of the punishment p. In some situations these can be given 
rough measures. To a first approximation one might set b, b,, ¢ and 
¢, proportional.to r and $ and y proportional to p. We may also in- 
troduce explicitly a variable D which is zero for the case of recogni- 
tion learning and one for recall. The values of w can be obtained for 
D=0Oand D=1. If these are w’ and w”, we may set w= (1 — D)w' 
+ Dw”, thereby combining both cases into a single expression. Simi- 
larly, instead of w) we may write Dw, so that it vanishes for D = 0 
and equals w, for D=1. The dependent variable c may be given by 
n—-wW—W,wand w, being the new variables. We thus have the two 
dependent variables w and w, determined as functions of the seven 
independent variables n, M,N,f,D,,r and p. 

For a preliminary investigation one might set 6 = 0, and intro- 
duce the same coefficient « in the expressions ¢(M), 7(M), and ¢,(M), 
although in ¢, and perhaps also in 7, one might set u = 0. It would 
also seem plausible to set b/f = b,/¢,. Thus if b = ar, Cf = ar, 
b, = ar, 0, = AszMgr/ Ay , p= Gap, n = Asp, U = A and h, = a; , we 
shall require seven parameters, @,, d) --- a; , to determine a curve in 
a nine dimensional space. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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BOOK REVIEW 


D’ARCY WENTWORTH THOMPSON. On Growth and Form. Second Edition, 1942. 
Ir + 1116 pp. Cambridge: At the University Press; New York: The Mac- 
millan Company. $12.50. 


The new edition follows its predecessor after 25 years and has grown from 
729 to 1097 pages. Two chapters, namely III and VII, have doubled in vol- 
ume, while others exhibit a correspondingly lower rate of relative growth. 

The purpose of the book is stated on page 15: “We want to see... how the 
form of living things . . . can be explained by physical considerations,” and in 
the Epilogue the author writes: “My task is finished if I have been able to show 
that a certain mathematical aspect of morphology is essential to (the) proper 
study of Growth and Form. ... For the harmony of the world is manifest in Form 
and Number.” The term “Form and Number” is used once more, and the Epi- 
logue ends with a eulogy of Henri Fabre who “saw in Number le comment et le 
pourquois des choses.” This change of “Growth and Form” into “Form and Num- 
ber” and finally into “Number” appears, to the reviewer at least, as significant. 
For the book does not deal so much with relation between Form and Growth, as 
with Growth and with Form. The bond between the two is frequently no more 
than a formal one, namely, the mathematical treatment of the subject. 

Thompson’s conception of form is of classical simplicity. “When we abstract 
. .. from the material to the form, then force is the appropriate term for... 
the causes by which these forms... are brought about.” Sir Oliver Lodge is 
cited: “You can move matter, it is the only thing that you can do to it.” 

The principles according to which Form is explained vary. Chapter V con- 
tains, apropos of surface tension and its influence on Form, an explanation of 
the method of maxima and minima and the principle of least action. Another 
viewpoint is developed in the discussion of the equiangular spiral: “The evidence 
of (a) constant relativity of Growth is of the essense of the equiangular spiral.” 
It is here that form is directly considered as a result of growth. 

Growth forms the subject matter of Chapter III, entitled “The Rate of 
Growth.” It occupies itself mostly with ontogenetical problems. In discussing 
brain weight and body weight, however, there appears a table illustrating that 
“the ratio of magnitude tends to change not only as the individual grows but also 
with change of bodily size from . .. one species to another,” and near the very 
end of Chapter III we find the statement that “changes of form are essentially 
the same . .. whether they be episodes in the life history of the individual, or 
mainfest themselves as the distinctive characteristics of the species of the race,” 
(the word race is here evidently intended to denote genus, not a subdivision of a 
species). “From one form to another there is but one road of transformation, 
and... it will proceed in the self-same way whether it occurs within the lifetime 
of an individual or during the long ancestral history of a race.” 

This shift from ontogenesis to phylogenesis is completed in the last, and prob- 
ably most famous chapter, “On Transformation.” The principle is well known: 


het 
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Of two (or several) forms to be compared, one is inscribed into rectangular Car- 
tesian coordinates. It is required to find the deformation (or deformations) these 
coordinates have to undergo in order to reproduce the other form (or forms). The 
chapter starts out with simple examples, but leads up to problems requiring cur- 
vilinear deformations. No attempt is made to handle these problems analytically. 
The diagrams are startling and fascinating, and it is all the more a matter of 
surprise that they have been comparatively little used by other authors. It still 
remains to be seen whether an improved technique could render the method more 
useful. The method obviously represents problems of conformal representations. 
Thompson cites a paper by Ikada and Kuwaori on conformal representation by 
means of elliptic integrals. A powerful tool could perhaps be forged out of the 
methods of these Japanese authors, which, written from the standpoint of a phy- 
sicist, are not immediately applicable to Thompson’s problem. 


Thompson is essentially a biologist who argues from the concrete instance to 
the general. This bent of mind is well illustrated by his discussion of cell division. 
It is interesting to compare Thompson’s treatment with that of Rashevsky, Lan- 
dahl and Young, which is given but scant justice in the book. Rashevsky starts 
out from the simplest assumption of the cell as a metabolizing system. Thomp- 
son watches under his microscope how a cell divides and gets immediately involved 
in karyokinetic spindles.. He does not get beyond a vague analogy with magnetic 
forces, culminating in the statement (in italics): “A body placed in the medium 
will tend to move towards regions of greater or lesser force according as its per- 
meability is greater or less than that of the surrounding medium.” 


The book contains a prodigious wealth of information., Written in a broad 
style and with frequent digressions from the main topic, it has, for these very 
reasons, a great and unique charm. On opening it at almost any place, one feels 
immediately that warmth which only civilized and generous persons are able to 
impart. And Thompson’s erudition goes far beyond the confines of biology. Ref- 
erences to the older literature are frequent, and philosophers and mathematicians 
are cited freely, be it in English, French, German or ancient Greek. 


The forms of animals and plants are universally explained on the basis of 
organic evolution. What is the author’s attitude toward this concept? On page 11 
we read: “Organic evolution has its physical analogue in the universal law that 
the world tends to pass from certain less probable to certain more probable states 

.. ” and a few lines further, this is described as the law of entropy. While 
Thompson states on page 14 that he is not concerned with the “weighty problem” 
of explaining “the long battle against the cold and darkness which is death,” he 
nevertheless appears to have opened the door to misunderstanding. 


Progress in organic evolution has been defined as increased independence of, 
and control over, the environment. This clearly makes evolution dependent of the 
environment which in its turn must be thought of as a function of entropy. In 
the light of the definition just cited, organic evolution would be rather dependent 
upon, than an analogue of, the law of entropy. To- reconcile entropy and or- 
ganic evolution is one of the most difficult problems of biology. Entropy is defined 
as probability and the law of entropy states that the final state of the world is 
probably predetermined. Does the analogy assumed by Thompson imply that, to 
cite but one concrete instance, homo sapiens is a more probable state of matter 
than a unicellular organism? There are, of course, many féwer men than there 
are microorganisms on this earth, therefore homo sapiens is less probable in the 
strictly statistical sense of that term. There is, as far as we can see, no difference 
in the probability of man and of amoeba in a thermodynamical sense. 
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It will readily be admitted, however, that the human organism is more com- 
plex than amphioxus or an amoeba. Complexity may be hard to define. The num- 
ber of motor patterns which an organism is capable of achieving may afford an 
objective measure of complexity. Thus fishes have but few motor patterns: loco- 
motion, steering, breathing, mastication, eye movements probably complete their 
inventory. Man, on the other hand, can move his head, his hands and his arms in- 
dependently from the rest of his body, can use his larynx for articulate speech 
and his facial musculature to express his emotions. Now if it is admitted that 
complexity is a less probable state for (living or dead) matter than simplicity, 
then it follows that organic evolution and the law of entropy are not analogues 
but opposites. “Life resembles the sailor who climbs up the mast of a sinking 
ship,” is the gist of a remark that the reviewer read years ago but failed to track 
down again. If organic evolution were a case of entropy, the battle against cold 
and darkness could never have been waged. 


In an addition to the second edition, the author cites Kant who states that 
the “axiom of casuality is denknotwendig, for experience itself depends upon its 
truth.” (p. 21). Elsewhere, (p. 268 et seq.) Thompson hints that “the great 
function of natural selection is not to originate but to remove,” and adds that 
“the world of things living .. . pursues its ceaseless course of creative evolution.” 
One is reminded of a passage of Whitehead’s “Adventures of Ideas” (p. 368/369) : 
“Decay, Transition, Loss, Displacement belong to the essence of the Creative 
Advance.” But Thompson is far from Whitehead, whose name does not even ap- 
pear in the index. As we just saw, Thompson follows Kant in his conception of 
the law of causality. What, then, is meant by creative evolution? Clearly, crea- 
tivity implies the negation of an unrestricted reign of causality, for it presup- 
poses spontaneity. It is interesting to see here how Thompson, though deeply 
steeped in mechanism, has still not completely avoided the terminology of an en- 
tirely incompatible position. 

Furthermore, the connection between ontogenetic and phylogenetic processes 
appears in practice to be far from “easy and legitimate”, as the author states. 
The problem of brain weight and body weight, adduced by the author, is a case 
in point. It is easy to find a comparatively simple functional relation between 
these two variables when the universe of discourse is defined, let us say, to be 
the class of adult mammals. But when we deal with the class of all develop- 
mental stages of a mammal, it is not possible to find a simple functional rela- 
tionship, since, as Scammon has shown, body-weight follows another type of 
growth than brain weight. Since this difficulty enters into most problems con- 
cerning mammals, any attempts to explain mammalian evolution by allometric 
growth are beset with difficulties. That many problems can be attacked in that 
way has been shown by Julian Huxley (Problems of Relative Growth, New York, 
1932) who dedicated his book to Thompson, 

The enthusiastic attempt by Huxley, Ray Lankaster, Gegenbaur and other 
fathers of modern morphology to understand the form of living things as a prod- 
uct of their phylogenesis has largely failed and been left unfinished. This may 
partially have been due to lack of material, but it is also due to insufficient analy- 
sis of the logical foundations of morphology. It would be preposterous to expect 
such an analysis in “Growth and Form”, but it is to be regretted that Thompson 
does not seem to be aware of its need and by unconsciously evading the issue may 
have discouraged progress in this field. 

In spite of its great charm (there is no other word for it), the book is not 
forceful. Partly, at any rate, this may be ascribed to Thompson’s failure to sub- 
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mit his problems to a logical analysis that could compare in rigor with mathe- 
matical analysis. The book is so full of stimulating thought, however, and con- 
tains such an abundance of examples that one can feel nothing but gratitude for 
Sir D’Arcy Thompson for having undertaken the considerable labor of revision 
and re-editing. Even if the reviewer disagrees with some of its tenets, he is well 
aware of the fact that these differences concern mainly questions about which 
there is not. and cannot be complete agreement. It is not a defect to have made 
statements open to attack, it is rather a merit to have stimulated thoughts on the 
fundamental problems of biology. 

GERHARDT VON BONIN 
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